This paper deals with the numerical treatment of singularly perturbed parabolic reaction-di, usion initial boundary value problems. Introducing a , tting parameter into the asymptotic solution and applying average , nite di, erence approximation, a , tted operator , nite di, erence method is developed for solving the problem. To ac-celerate the rate of convergence of the method, Richardson extrapolation technique is applied. The consistency and stability of the proposed method have been established very well to ensure the convergence of the method. Numerical experimentation is carried out on some model problems and both the results are presented in tables and graphs. The numerical results are compared with , ndings of some methods existing in the literature and found to be more accurate. Generally, the formulated method is consistent, stable, and more accurate than some methods existing in the literature for solving singularly perturbed parabolic reaction-di, usion initial boundary value problems.